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Abstract 

We are interested in the experimental characterization of the Hofstadter butter- 
fly by means of acoustical waves. The transmission of an acoustic pulse through an 
array of 60 variable and resonant scatterers periodically distribued along a waveg- 
uide is studied. An arbitrary scattering arrangement is realized by using the vari- 
able length of each resonator cavity. For a periodic modulation, the structures of 
forbidden bands of the transmission reproduce the Hofstadter butterfly. We com- 
pare experimental, analytical, and computational realizations of the Hofstadter but- 
terfly and we show the influence of the resonances of the scatterers on the structure 
of the butterfly. 

PACS. 43.20.Fn-Scattering of acoustic waves 

PACS. 43. 20. Mv- Waveguide, wave propagation in tubes and duct 

PACS. 43.58.Gn-Acoustic impulse analyzers and measurements 

1 Introduction 

The Hofstadter butterfly is well known since the work of Hofstadter on the electron 
transmission through a two-dimensional ordered lattice perturbed by a perpendicularly 
applied uniform magnetic field [JJ. Few experimental attempts have been made to 
observe signatures of this phenomenon in two-dimensional electron systems 121 [3] 14]. 
Recently, an experimental electromagnetic realization of the Hofstadter butterfly was 
done by studying the transmission of microwaves through an array of 100 scatterers 
inserted into a waveguide when the modulation length of this unidimensional periodic 
lattice changed [5]. These authors used cylindrical scatterers introduced into a rect- 
angular microwave waveguide. In this study, the stopbands were produced by Braggs 
reflection on the transmission spectra and the illustration of the Hofstadter butterfly 
exhibited self-similar structures. 
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In the present paper, we study the transmission of an acoustic pulse through a lat- 
tice composed of an assembly of variable Helmholtz resonators that are the scatterers 
||6l. A particularity of this work is the use of scatterers that can exhibit resonant scatter- 
ing at wavelength much larger then the scatterers (well known examples of Helmholtz 
resonators are the bottle of wine that have resonances for frequencies corresponding 
to wavelengths of the order of the meter). The resonators are uniformly distributed 
along a cylindrical waveguide. Each Helmholtz resonator is constituted by two cylin- 
drical tubes which play the role of neck and cavity of the resonator The length of the 
Helmholtz cavity is a variable parameter and the lattice spacing constant d between 
each resonator is fixed. 

The propagation of sound wave into a tunnel with a periodic array of Helmholtz 
resonators has been already studied Q and different kinds of stopband appear in the 
transmission coefficient. In addition to Bragg stopbands caused by the spatial periodic- 
ity of the lattice, resonances of the Helmholtz resonator create other "stopping bands" 
called scatterer stopbands which are non-existent in |5|. We show that a realization of 
the Hofstadter butterfly is obtained experimentally with acoustic waves. The effects of 
the scatterer stopbands are examined. Experimental work is completed by an analytical 
and a numerical analysis, and all show self-similar structure. 



2 Theory of the propagation into a lattice 

We study a monochromatic pressure wave p{x,t) with the form p{x,t) — p{x)e^^ 
where co is the wave pulsation. A cylindrical waveguide with Helmholtz resonators 
connected axially constitutes the one-dimensional lattice. In this lattice, the pressure 
p{x) is solution of the following equation flSl, 

^-^^+k^p{x)=Y,^{x-Xn)Onp(,x) (1) 

where k = (x)/c, a„ = —jOip^Yn with Si and S the derivation and the main tube areas 
respectively (fig. [T]i, p is the density of air and c is the celerity of the wave. F„ is the 
admittance of the n''' resonator placed atx = x„, defined by Y„ = v{xn) / p{xn) = v„/ pn 
where p{xn) and v(jc„) are respectively the pressure and the acoustic velocity ?Ax — x„. 
The solutions of eq. ([T]) are found by the transfer matrix method. In the region x„ < 
X < Xn+i, which determines the n'^' cell, the solutions of eq. ([T]) are separated into two 
plane waves with opposite propagation direction associated with the amplitudes A„ and 
B„: p{x) =A„e^''^(^--'-«)+B„e-^''^(-»'--'^"). 

The propagation through one cell is described by 

^"^'^-mJ^"\ = \ - X':m. ..^%,..-m. 1(7 ) (2) 




with dn — Xn+i —x,,. For a finite lattice which is composed of resonators, the wave 
propagation is described by {A„B„y = Yl^=]_M„{AoBoy , where {AqBo)' are the ampli- 
tudes of the inlet of the lattice. 
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Figure 1 : Description of the Helmholtz resonator. 



The resolution of the eigenvalue problem for the matrix M„ allows us to find the 
structure of the transmission (allowed and forbidden bands). By setting A,* — e^^^"^^ for 
the eigenvalues, the dispersion relation takes the form 

cos{qd) = cos{kd) + ^ sm{kd), (3) 

where q is called the Bloch wave number The transmission of the lattice is deter- 
mined by the value of cos(qd) : if cos{qd) belongs to [—1, 1], the wave propagates 
into the system; in contrast, when cos{qd) belongs to ] — — 1] or [1 , +°°[, the wave is 
exponentially damped and the lattice is completely opaque to the wave. 



3 Determination of Helmholtz resonator admittance Y 

To determine the impedance of a Helmholtz resonator, we need to a relation pressure 
p and velocity v between the entrance (E) and the end (F) of the resonator (fig. [1). 
Applying classical lossless plane wave technique it is found that 



PE \ _ ( cos{kli) jZi.sm{kli) \ / 1 jzikAl 



VE J \ ^sin(fc/,) cos{kli) j y 1 
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cos{klc) jZ^.sm{klc) 
4^sin(fcZc) cos( klc) 




) 



(4) 



where ZJ. and Z^ are respectively the characteristic impedance of neck and cavity, Z, 
and Ic are respectively the length of the neck and the length of the volume cavity, and 
subscripts E and F refer to entrance and end of the resonator. Al represents the added 
length to the neck due to the discontinuity of the sections. This length takes generally 
the form f9l A/ A/,- + Al,- = '"Kt + |e + where n is the radius of the neck. 

Because the Helmholtz resonator is closed, the boundary condition at the end of the 
resonator is vp = 0. Eliminating pf in the eq. (01, the acoustic impedance Z = pe/ve 
seen from the entry (at the point E) is given by the following relation : 



Because of resonances, it is very important to take into account the losses. The 
theory with losses OO] is used in the eq. ^ with k = ^(1 + 7(1 + (y— l)x) ™d 
Z^. = pc(l + ^(1 — (y— l)x) by setting s = r/5 where 5 is the viscous boundary layer, 
X = \/P^- with Pr the Prandtl number, p = (1 — /)/ \/2 and where r is the radius of the 
tube considered. 

4 Description of the experimental apparatus 

In the experiment, a cylindrical waveguide with an inner radius r = 2.5 10^^ m was 
used. 60 Helmholtz resonators were periodically located along this waveguide and 
the lattice spacing between each resonator was d — 0.\ m. Each resonator was com- 
posed by a neck (cylindrical tube with an inner area Si = 7,85 .10^^ and a length 
/,■ = 2.10-2 m) and a variable length cavity (cylindrical tube with an inner area Sc = 
7,85.10-^7i2 and a maximum length ~ 16.5 .10 ^ m) as described in the fig. [T| The 
sound source was a loudspeaker with special design placed at one end of the main tube 
. Two microphones (BK 4136) measured the pressure at each end of the cylindrical 
waveguide. We measured in the frequency range where only the first mode can prop- 
agate, starting from Hz up to the first cutoff frequency /eoi =4,061 kHz, where the 
propagation of the first higher mode becomes possible. 

For the experimental realization of the butterfly a periodic modulation of the lengths 
of the cavity resonators was applied with the period length as a parameter. The modu- 
lation of this lattice is given by the variation of the cavity length of the n''' resonator 



where a = ^, with m e N. We replace this cosinusoidal variation by a rectangular 
one (in the same manner as ||5l): 



Z 




(5) 



y cos{klc) sm{kli) — ^ sin(fcZ,) sin{klc) + cos,{kli)sin{klc) 



l„ = /tCos(27tna — ao) 



(6) 



{ 



for cos(27ina-ao) < 
Ic for cos(27tna— ao) > 0. 



(7) 
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In the experiment we choose Ic — 16.5 lO^^m and ao = 0. It is important for the 
experiment to use a rectangular setup and not a sinusoidal otherwise the stopbands 
would destroy the butterfly. 



To investigate the experimental results, we used a time-frequency method. It maps the 
time domain signal (the square of the amplitude) in the time-frequency plane. In this 
way, it is possible to know the arrival time of each frequency present in the input signal 
and to determine a "transfer function" of the acoustical system. The Wigner-Ville ifTTl 
distribution defined as 



where Zx{t) is the analytical signal lfT2l (associated to the real signal z(f )) is certainly 
the most widely studied time-frequency methods. To reduce the cross terms (interfer- 
ences) between different components of the signal, we used the Pseudo Wigner-Ville 
distribution PW given by : 



where h{z) is the temporal rectangular window with a length x. This method gives the 
acoustical response of the lattice in the time-frequency plane by using a source pulse. 
We determined the transmission spectra by detecting the maximum of energy for each 
frequency upstream and downstream the lattice. The transmission was deduced from 
the "transfer function" (in energy) defined as the fraction of transmitted and incident 
energy. 

Using a pulse source, this method avoids seeing, in the transmission spectra, os- 
cillations phenomena caused by the finite number of cells 1 13 1. The pulse source pre- 
vents this phenomena, which would destroy the Hofstadter butterfly structure, when the 
propagation time along a cell is greater than the pulse duration 1 14|. The fig. |2]shows 
a comparison between these two sorts of processing method : a classical way (Fourier 
transform) using a shirp source and the previous one using the pseudo distribution of 
Wigner-Ville and a pulse source. It illustrates the oscillations in the transmission coef- 
ficient and shows the influence of the source on the transmission spectra. 

6 Analytic determination of the Hofstadter butterfly 

The periodic rectangular modulation implies only two values for the cavity length of 
each Helmholtz resonator. When this length vanishes, the resonator effects on the 
propagation are completely invisible in the frequency range of interest. When a is 
written as a = 1/m, the lattice is periodic with a period length d' = d/a = md. In this 
case, analytical points of the Hofstadter butterfly can be found in the (/, a) plane. 



5 Signal analysis 
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Figure 2: Comparison between two different signals processing for a — 0.0125. (a) 
Use of a Fourier transform and a shirp source . (b) Use of Wigner-Ville distribution 
and a pulse source. 



This new period d! is used in the dispersion relation (eq. (O) which becomes 

cos(/tt/') + ^sin(;tt/') =±1- 

The solutions of this equation gives some points of the line which is the boundary 
between the passbands and the stopbands (the stopband are determined when the dis- 
persion relation (O is greater than 1 or lower than —1). By setting tan(9) = — the 
dispersion relation is written as cos(^^ + 9) = ±cos(9) and the solutions are 



d d 
k— = nn — 29 and k— — rni with n e N. 
a a 



(8) 



The effects of the scatterers are present in the term 9 and the resonance frequencies 
imply its divergence. These divergences perturb branches structure of the butterfly 
around these resonance frequencies and "break" the butterfly structure as presented in 
the fig. |3]and|4]for / = 283 Hz, / = 1 120 Hz, / = 2240 Hz and / = 3360 Hz (each 
band corresponds to a resonance of the resonator). The resolution of the eq. (O gives 
different points on the (/, a) plane for a's values corresponding to p = 1 and m = 2 up 
to 20 which are reported on the fig. [3j). 
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Figure 3: (a) Experimental transmission spectra for a periodic arrangement of scatterers 
with a ranging from to 1 with 50 values between and 0.5. The upper part is 
obtained by reflection, (b) Zoom of the experimental transmission spectra for a periodic 
arrangement of scatterers with a ranging from to 1 between [0 : 2000] Hz. The white 
points show the analytic results from the eq. (O. 

7 Results 

The transfer function between the end and the beginning of the lattice, in the frequency 
range [0,4000] Hz which is corresponding to two 7t-Bragg bands, was measured. We 
used 50 different values for the parameter a between and 0.5. 

The fig. |5] presents the transmission in the periodic case (all the resonators are 
identical) and for three different values of modulation length. The periodic case (fig. 
|5h) is obtained for a = and shows two sorts of forbidden bands in the spectra : the 
first and second are due to the scatterers (Helmholtz resonators) for 283 Hz and 1120 
Hz and the third is due to the spatial periodicity (7t-Bragg stopband) at 1700 Hz. 

The periodic modulation of the lattice involves the arrival of forbidden subbands in 
the transmission spectra (fig. |5}3, c and d) which split the transmission band (allowed 
band) into subbands (marked with a arrow in the figure). The number of these subbands 
increases with m. The new period of the lattice is d' ~ md so the frequency location of 
the subbands are now ft—c/ {2md) : 1 subband appears for p — 1 and m = 2 (fig. |5j)) at 
850 Hz and theoretically 10 subbands for p = 1 and m = 10 (fig. |5}J) every 170 Hz (this 
argument is not satisfactory for /:> 7^ 1 as in the fig. |5}; where p ~ 13 and m — 40). The 
attenuation phenomenon added to the finite dimension of the lattice and the presence 
of stopbands due to resonance's scatterers induce that some subbands are missing in 
the spectra presented in the fig.H] 

To obtain a view of the Hofstadter butterfly, the transfer functions for different 
values of the modulation parameter a are plotted together on a plane building with 
the frequency (abscise axe) and a belonging to [0, 1] (ordinate axe). The spectra are 
converted to a grey scale (fig. |3]and|4]i and the part between [0.5, 1] is obtained by 
reflection. 

In fig. |3^ two 7t-Bragg bands are seen and two Hofstadter butterflies are complete. 
The structure of the butterflies is broken by the resonances of the resonators and some 
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Figure 4: Numerical result for the transmission spectra for a periodic arrangement of 
scatterers with a ranging from to 1 with 100 values between and 0.5. The upper 
part is obtained by reflection. 

"oscillations" due to the signal processing are observed. The fig. ^ presents a zoom 
of the fig. [3^ for a [0 : 2000] Hz frequency range. In this figure, the structure of the 
Hofstadter butterfly are identified and can be compared to the simulation (fig. |4|i ob- 
tained with 100 values of a. Self similar structures are observed and a lot of subbands 
are perceptible. The presence of scatterers stopbands seems to stop the effects of the 
lattice modulation in the transmission : for example, the main subband is just shifted 
during the both stopbands caused by Helmholtz resonators and the butterfly structure 
is affected by this phenomenon. The white points on the fig. [Sj) present the analytical 
results calculated with the eq. dHJ. They are in agreement with the experimental results 
and they predict sufficiently the effect of the scatterer resonances. 

The self-similar structure of the butterfly is also observed. The simulated butterfly 
(fig. nil, building with 100 different arrangements of the scatterers, exhibits more self- 
similar structures than the experimental results. The depth of the self-similar structures 
is depending on the number of modulation used to map the Hofstadter butterfly. The 
simulation allows to observe with a high definition some fine details of the Hofstadter 
butterfly and to analyze more precisely the effects of the scatterer resonances. It seems 
that the presence of scatterer stopbands (the first and the second stopbands) breaks the 
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Figure 5: Transmission function for four different values of modulation lengths, (a) 
a = 0. (b) a = 0.5. (c) a = 0.325. (d) a = 0.1. Some subbands due to the modulations 
are marked by arrows. 

modulation influence by diverting the different branches of the butterfly. It implies that 
the structure is self-similar only in the region between the scatterer stopbands. This 
phenomena is also observed with the experimental results for the main branch of the 
butterfly (fig-Wp)- 
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